Giant Spin Motive Force due to Rashba Spin-Orbit Coupling 
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Magnetization dynamics in a ferromagnet can induce a spin-dependent electric field through spin 
motive force. Spin current generated by the spin-dependent electric field can in turn modify the 
magnetization dynamics through spin-transfer torque. Whereas this feedback effect is negligibly 
weak in conventional situations, we report that the spin-dependent electric field in strong Rashba 
spin-orbit coupling systems can be orders of magnitude stronger and generate spin current large 
enough to modify the magnetization dynamics significantly. This effect should be relevant for 
device applications based on ultrathin magnetic layers with strong Rashba spin-orbit coupling. 



Similar to the mutual induction between electric and 
magnetic fields through the Faraday and Maxwell's laws, 
spin current and magnetization induce the dynamics of 
each other in magnetic systems through spin-transfer 
torque (STT) [Hlf and spin motive force (SMF) j^Q. 
Through many years of extensive study 0, 0, H| , it has 
been demonstrated that STT is a powerful tool to induce 
the magnetization dynamics in ferromagnetic nanostruc- 
tures. In contrast, SMF has received much less atten- 
tion since SMF is very weak and thus inefficient to 
generate the spin current. 

In this Letter, we demonstrate that SMF can be orders 
of magnitude enhanced in systems with strong Rashba 
spin-orbit coupling (RSOC). RSOC arises generically 
when structural inversion symmetry is broken [l2| . Re- 
cent experiments on ultrathin (~ 1 nm) ferromagnetic 
layers with the stron g st ructural inversion asymmetry 
(such as Pt/Co/AlO,, [I3,[i3 and Ta/CoFeB/MgO [H) 
observed large effective magnetic fields predicted by 
RSOC theories 0, [ljj ■ A thin ferromagnetic layer in 
contact with topological insulators may also have 

strong RSOC. In such magnetic systems, magnetization 
dynamics can induce large spin current through SMF, 
and this spin current can in turn modify the magne- 
tization dynamics through STT. Thus even for purely 
magnetic- field- driven magnetization dynamics, STT can 
have sizable magnitude because of large spin current 
generated through SMF. We also propose a SMF-based 
method to measure the strength of RSOC in magnetic 
systems. This method allows for unambiguous distinc- 
tion between RSOC effects and other effects [20j . which 
may be difficult to distinguish through other methods. 
Thus SMF can also be a useful tool to quantify RSOC. 

RSOC can be described by the Rashba Hamiltonian 
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where the vectors cr and p are the Pauli matrix and 
the momentum, respectively, z is the unit vector along 
the inversion symmetry breaking direction (perpendicu- 
lar to ferromagnetic layer), and cxr is the Rashba con- 
stant. The total Hamiltonian H of a conduction electron 



then becomes Hq + Hr, where Hq = p 2 /2m e — J ex cr ■ m. 
Here m e is effective mass of conduction electrons, J ex 
(< 0) is the exchange coupling energy, and m is the 
unit vector along the local magnetization, which is in 
general time- and position-dependent. An insight into 
the RSOC effect on SMF can be gained from the ve- 
locity operator v = [r, H]/ih = p/m e + v an , where the 
anomalous velocity v an = [r, HrI/Ht, = o.rz x cr jh arises 
due to RSOC. When the exchange energy is sufficiently 
larger than RSOC, cr for majority (minority) electrons 
will be almost anti-parallel (parallel) to m. It is then 
evident that the magnetization dynamics dm/dt induces 
the acceleration dv an /dt and the effective electric field 
— (m e /e)dv an /dt. This heuristic calculation results in the 
RSOC contribution to the spin-dependent electric field 
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where ± applies to the majority and minority electrons, 
respectively. The total spin-dependent electric field E' ± 
then becomes E± + E 
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where 
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represents the conventional spin-dependent electric field 
examined in previous literatures [7|, |21| in the absence 
of RSOC. Supplementary information [22| provides more 
rigorous derivations of Eq. ^ based on the SU(2) gauge 
field formalism Q and the spin misalignment force for- 
malism [23j. 

We estimate the relative magnitude of E^ soc with re- 
spect to E±; E^ soc is of the order of (aRm e /eh)uj, where 
uj is the characteristic angular frequency of the magneti- 
zation dynamics, and E± is of the order of (H/e)(uj/L), 
where L is the characteristic length of a magnetic struc- 
ture [such as domain wall (DW) width] producing E±. 
The relative ratio becomes aRm e L/h 2 ~ cxrL x 10 38 
(kg/J 2 - s 2 ) for m e — 9.1 x 10~ 31 kg. For the value 
a R = 1CT 10 eV-m reported for Pt/Co(0.6 wn)/MO x [H 
and for L — 20 nm, this ratio becomes about 30. Thus 
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FIG. 1: (color online) Voltage produced by a precessing DW 
in a nanowire. For simplicity, the DW motion along the 
nanowire is assumed to be suppressed by a notch in the 
nanowire. (a) Schematic structure of DW. DW is assumed 
to precess and change its structure periodically between the 
Bloch DW (tp = or tt) and the Neel DW (ip = tt/2 or 3tt/2) 
with uj = dip/dt = 27rxl00 MHz. (b) Voltage between the 
two ends of the nanowire due to E± . This voltage does not 
depend on time t and its magnitude is comparable to the 
value reported in Ref. 0. (c) Voltage due to E± soc , which 
oscillates with t. Here an — 10 -10 eV-m [l3( |. P = 1, and the 
DW width of 20 nm are assumed. 



EjRSOC can be more than an order of magnitude larger 
than E±. 

For magnetic materials with nonzero spin polariza- 
tion P, spin-dependent electric fields can be measured 
through conventional electric voltage measurement [Toj . 
Figure [T] compares electrical voltages produced by a pre- 
cessing DW in a nanowire. The voltage due to E^ soc 
[Fig. [He)] exhibits time-dependence different from that 
due to E± [Fig. Ufb)] because of the functional form dif- 
ference between Eqs. ([2]) and <j3j> - Note that the oscilla- 
tion amplitude in Fig. [He) is about 80 times larger than 
the dc value in Fig. D^b). 

In regard to the value of «_r, 10 -10 eV-m for Pt/Co(0.6 
nm)/A102; [l3] is not exceptional. The an value in the 
range (0.4 — 3) x 10 -10 eV-m was reported by photo- 
electron spectroscopy measurements 24|, |25[ for a thin 
nonmagnetic metal layer in contact with a heavy atomic 
element layer. Since magnetism does not suppress an, 
its value can be in a similar range for proper com- 
binations of a thin magnetic layer in contact with a 
heavy atomic element layer. A recent effective field mea- 
surement result 15[ for Ta/CoFeB(f nm)/MgO implies 



aji 0.2 x 10~ 10 eV-m, which falls close to the com- 
mon range. For this particular multilayer, it was sug- 
gested [15J that the contact with the oxide layer MgO 
may also be an important source of olr. Still another 
candidate system is a thin magnetic layer in contact with 
topological insulators. For instance, BiTel has a large 
bulk RSOC with a R = 3.8 x lO" 10 eV-m [H] and was 
predicted [l9| to become a topological insulator under 
pressure. 

Large spin-dependent electric field implies large spin 
current. The spin current density generated by SMF be- 
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<t c (E+ + E+ uo ) [26|, where is the longitudinal elec- 
trical conductivity of majority (minority) electrons and 
Cc = o"t + a l ■ Since E^ soc is much larger than E + and 
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a magnetic layer. Then for a typical in-plane charge con- 
ductivity <t c = 10 7 ■ m -1 of metallic ferromagnetic 
layers and a R = (0.2 - 3) x 10~ 10 eV-m, J^ MF has mag- 
nitude (1.8 - 27)w A-s/m 2 . For u = 2ir x 5 GHz, which 
is a typical angular frequency scale for the magnetiza- 
tion dynamics [8j in magnetic tunnel junctions, we find 
|jS MF | = (0.6 - 8.5) x 10 11 A/m 2 . Thus RSOC allows 
SMF to generate huge spin current density from typical 
magnetization dynamics. 

Large jf MF generated by magnetization dynamics im- 
plies that magnetization dynamics itself can be signifi- 
cantly modified by j|? MF through STT. Next we exam- 
ine this feedback effect. In conventional situations where 
SMF is negligible, the magnetization dynamics is de- 
scribed by the Landau-Lifshitz-Gilbert (LLG) equation, 
dm/dt — 7m x H c ff + agm x dm/dt + T(J S ), where T 
represents the STT and depends on externally supplied 
spin current density J s . Here 7 is the gyromagnetic ra- 
tio, ac is the Gilbert damping parameter, and H c g is the 
sum of an external magnetic field and effective magnetic 
fields due to magnetic anisotropy and magnetic exchange 
energy. To examine the feedback effect, one simply needs 
to replace T(J S ) by T(J S + jf MF ) = T(J S ) + T(J™ F ) to 
obtain the modified LLG equation, 

^ = -7mxH cff +a G mx^+T(J s )+T(J s SMF ). (4) 

Note that STT T now has two spin current sources, J s 
andJ s SMF . Thus even when there is no externally sup- 
plied spin current (J s = 0), STT still affects the magne- 
tization dynamics as long as Jg MF is not zero. Therefore 
STT becomes relevant even for purely field-driven magne- 
tization dynamics. Moreover STT can become relevant 
even when both external field and external current are 
absent if m is initially far from its stable configuration. 

To gain an insight into roles of the feedback STT 
T(Jg MF ), we express it in the following form, 



T(jf* F ) = m x V 



dm 



(5) 
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This form is natural since T is orthogonal to m [thus 
T(J S SMF ) = mx (function of J S SMF )] and J S SMF is propor- 
tional to dm/dt [thus (function of J™ F ) = V ■ dm/dt}. 
Here T> will depend on m and be a 3 x 3 matrix in gen- 
eral. Although V is not a constant, the structural sim- 
ilarity between Eq. ([5]) and the Gilbert damping torque 
acrn x dm/dt implies that T> may be interpreted as a 
correction to the Gilbert damping parameter ag. Thus 
roles of the feedback STT can be analyzed by using this 
damping correction picture. 

After explicit calculation, we find that the matrix ele- 
ments of the 3x3 matrix V are given by 



Vij — V*^2 (A fe j + a R e 3ki ) (X kj + aR€ 3 kj) 
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where X ki = (m x dm/dx k ) v r\ = fj, B Ha c /2e 2 M s , a R = 
2a R m e /H 2 , /Jb(> 0) is the Bohr magneton, M s is the 
saturation magnetization, and €ij k is the Levi-Civita per- 
mutation symbol (ei23 = 1)- To derive Eq. ([6]), we 
simply combine the spin-dependent electric field equa- 
tions [Eqs. ((2J and ([3])] with known contributions to 
STT. To be specific, the adiabatic STT contribution 
[oc (J s • V)m] and the field-like STT contribution [oc 



a R m x (z x J s )] [16l . |17I | are taken into account in the 



calculation whereas the nonadiabatic STT contribution 
and recently discovered Slonczewski-like STT contribu- 
tion [27l - l29l | are ignored since the latter two contributions 
are smaller than the former two. 

The importance of the feedback STT can be estimated 
from the relative magnitude of T> with respect to etc- 
Here qq includes all contributions other than the SMF 
contribution. The intrinsic bulk Gilbert damping pa- 
rameter is of the order of 0.01 for typical metallic fer- 
romagnets. In a thin magnetic layer, this intrinsic value 
is enhanced by the conventional spin pumping mecha- 
nism [3(1 HH and the enhanced value is estimated to 
0.1/cZf (nm), which is 0.1 for thin magnetic layer of 
thickness dp = 1 nm. Thus for the feedback effect to 
be a relevant factor for the magnetization dynamics of a 
thin magnetic layer, T> should be comparable to or larger 
than 0.1. When RSOC is absent (a R = 0), V reduces to 
the result in Ref. 2l| and is of the order of rj/L 2 . For 
cr c = 10 7 rr 1 • m" 1 , M s = 10 6 A/m and L = 20 nm, 17 
is 0.2 nm 2 and rj/L 2 is 0.0005. Thus except for special 
situations with L < 1 nm, T> becomes much smaller than 
0.1 and the feedback effect is negligible. 

When RSOC is strong, on the other hand, the largest 
contribution to T> arises from the term na 2 B y \ en k.jtzk-i 
and thus V becomes of the order of a' G = r/a 2 R [32j . For 
a R = 10" 10 eV- m [H, & R is 2.6 nm" 1 and a' G is 1.4. 
Since a' G ^> oiq, one of evident effects of the feedback 
effect is to enhance the effective damping significantly. 
The second largest contribution to T> comes from the 
cross-term r]a R J2k( x ki£3kj + X kj e 3ki ). It is of the order 
of r\a R / L, which is 0.03 for L = 20 nm. This contribu- 
tion is somewhat smaller than aa but since it is inversely 



proportional to L, it also becomes relevant for localized 
magnetic structures with L < 5 nm. Thus another ef- 
fect of the feedback effect is to introduce the length-scale 
dependence to the effective damping. 

For deeper understanding of the feedback effect, how- 
ever, one should go beyond this analogy to the sim- 
ple form of the Gilbert damping torque agin x dm/dt 
and address the detailed structure of the damping ma- 
trix T> in Eq. ([BJ. Extensive discussion on implications 
of the damping matrix structure will be presented else- 
where. Here we present one simple implication contained 
in the largest contribution a' G J2 k e 3fei £ 3fej j which reduces 
to a' G {5ij— S^iSsj) after simple algebra. This contribution 
makes the effective damping anisotropic since Ylji^ij ~ 
83i53j)(dm/dt)j becomes (dm/dt)i when dm/dt is per- 
pendicular to i-direction, and vanishes when dm/dt is 
parallel to z-direction. For a magnetic layer with the per- 
pendicular magnetic anisotropy, this anisotropy can be 
tested in experiments by measuring the effective damp- 
ing in two different ways. When the effective damping 
is measured by ferromagnetic resonance, dm/dt remains 
essentially perpendicular to z-direction and thus the ef- 
fective damping becomes ag + a' G . The effective damp- 
ing may be measured alternatively by using the rela- 
tion that the field-driven DW velocity below the Walker- 
breakdown threshold is inversely proportional to the ef- 
fective damping. This measurement should produce a 
smaller effective damping value than the ferromagnetic 
resonance since dm/dt is along the indirection at the 
DW center. From a simple calculation based on the 
Thiele's collective coordinate description [33| of the DW 
configuration in terms of the DW position and DW tilt- 
ing angle, we find that the effective damping becomes 
a G + a' G /i- In this calculation, we assume that the DW 
width is sufficiently large so that the second largest con- 
tribution to V may be ignored. Thus in strong RSOC 
systems, where a' G 3> olq, we predict that the effective 
damping for the ferromagnetic resonance is about fac- 
tor 3 larger than that for the DW motion. By the way, 
Ref. [l3[ determined the effective damping for Pt/Co(0.6 
nm) / AlOz from the field-driven DW motion measure- 
ment and found 0.5, which is in reasonable agreement 
with our prediction ag + a' G /3 ~ «g/3 ~ 0.5 for the 
system. 

Lastly we propose an experimental setup to test the 
RSOC-enhanced SMF directly. As shown schematically 
in Fig. [21 when a uniformly magnetized single mag- 
netic domain precesses due to external magnetic field 
B{t) = B shout, E^ soc will generate electric voltages V x 
and V y . Since E^ soc is position-independent [Eq. ©], 
these voltages should be proportional to the spacing d 
between electrodes. Also since the direction of E^: soc ro- 
tates within xy-plane as dm/dt rotates within xy-plane 
(due to the perpendicular magnetic anisotropy), the volt- 
ages will be oscillatory with the oscillation phases of V x 
and V y different from each other by about 90°. These fea- 
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FIG. 2: (color online) Measurement scheme of the RSOC ef- 
fect on the SMF for a uniformly magnetized layer with the 
perpendicular magnetic anisotropy. 



tures are qualitatively different from other effects that 
may affect this measurement. The conventional spin 
pumping [3(| may also contribute to the electric voltages 
if the two electrical contacts with the electrodes are not 
identical, as demonstrated in Ref. (34|. These contribu- 
tion is however independent of d and thus distinguishable 
from the RSOC-enhanced SMF contribution. Spin Hall 
effect in contacts or the neighboring heavy metal layer 
can also generate a contribution. This spin Hall con- 
tribution is however different from the RSOC-enhanced 
SMF contribution since z x m determines the direction of 
the voltage gradient in case of the spin Hall contribution 
whereas z x dm/dt does in case of the RSOC-enhanced 
SMF contribution. Recalling that m and dm/dt are or- 
thogonal to each other, these two contributions differ by 
the oscillation phase of 90°. Recently it was reported [20( 
that some experimental results [35| , which was previously 
attributed to RSOC, can arise instead from spin Hall ef- 
fect in a neighboring heavy metal layer in contact with 
a thin magnetic layer. The electrical voltage measure- 
ment in Fig.[2]can be a very decisive experiment to verify 
the RSOC-enhanced SMF and also RSOC itself free from 
such controversies. The Rashba constant an can be de- 
termined from the proportionality constant between the 
voltages and d. 

In conclusion, we demonstrated that SMF can be or- 
ders of magnitudes enhanced in strong RSOC systems. 
The enhanced SMF is strong enough to modify the mag- 
netization dynamics significantly. RSOC-enhanced SMF 
may affect performance of various device applications 
based on ultrathin magnetic layers 35|-37| since RSOC 
effects tend to get larger in thinner systems. Thus SMF 
is not a weak effect of purely scientific concern any more 
but instead a technologically relevant effect that should 
be taken into account in future studies. As a final remark, 
during the revision of this manuscript, we were informed 
that two other groups [HEl] also obtained Eq. 
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